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PATENT 

Attorney Docket No. 6514 
ROD PUMP CONTROL SYSTEM INCLUDING PARAMETER ESTIMATOR 



CROSS REFERENCE TO RELATED APPLICATIONS 
[0001] This application claims priority of 

provisional application serial number 60/414,197, 

10 entitled "Rod Pump Control System Including Parameter 
Esitmator" , which was filed on September 27, 2 002, and 
provisional application serial number 60/429,158, 
entitled " Sensorless Control Sustem For Progressive 
Cavity and Electric Submersible Pumps", which was 

15 filed on November 26, 2002, and is related to 

application serial number entitled "Control System For 
Progressing Cavity Pump", which was filed on September 
5, 2003, and application serial number entitled 
"Control System For Centrifugal Pumps", which was 

20 filed on September 5, 2003, which was filed on 

September 5, 2003, which four patent applications are 
hereby incorporated herein by reference. 

BACKGROUND OF THE INVENTION 
25 [0002] Field of the Invention -- The present 

invention relates generally to control of rod pumps 
for oil and gas wells, and in particular to methods 
for optimizing the operation of rod pumps using 
parameter estimation . 
30 [0003] The load upon and position of the rods that 

drive downhole pumps are important parameters for 
control, monitoring, and protection of the artificial 
lift system used in oil and gas production. Existing 
methods of measuring these parameters involve the 
35 mounting and use of external instruments such as 
strain gauges, load cells, and position transducers. 
The need for these additional devices increases the 
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cost and complexity of the pumping system and reduces 
system reliability. Generally, AC induction motors 
drive rod pumping systems . 

[0004] One method for determining rod load or force 

5 is disclosed in United States Patent No. 4,490,094 

(the '094 Patent). With this method, motor velocity 
is determined during a complete or predetermined 
portion of a reciprocation cycle and the results are 
used to compute one or more parameters of pumping unit 
1 0 performance . 

[0005] However, determination of rod load PRL± on 
an ith revolution of the prime mover rotor depends on 
knowing the position of crank for computation of a 
torque factor TFi according to the equation (1) : 

15 

PRLi= n.Ti+m.sin(q> i +P) -RIT+A1T (1) 

TFi 

[0006] Because the torque factor TFi appears in the 

denominator of the equation, special care must be 

2 0 taken in deriving the torque factor Tf± and in using 

it in the computation to avoid dividing by zero or by 
small numbers that would distort the result. 
Moreover, the '094 Patent does not disclose how to 
estimate crank position. 
25 [0007] United States Patent No. 5,252,031 (the '031 

Patent) discloses a method for monitoring a rod pumped 
well to detect various problems. The method uses 
measurements made at the surface to calculate a 
downhole pump dynamometer card. This downhole pump 

3 0 dynamometer card is useful in detecting various pump 

problems and controlling the pumping unit. The method 
involves finding rod position from motor revolutions, 
a reference switch and pump geometry. This method 
requires setting up look-up tables. 
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[0008] In addition, the methods disclosed in both 

the '031 Patent and the '094 Patent employ a sensor to 
detect a rotation of the motor shaft. Because of the 
ratio between motor and pump rotations, this method 
5 can produce numerous sample points per stroke of the 
pump. However, the time between motor revolutions to 
get motor velocity as well as sample other parameters, 
such as motor current, is a function of pump speed and 
is not suitable for precise monitoring of the pump 

10 operation. In addition, the method of determining 
motor torque relies on a look-up table of steady-state 
motor operation rather than a true dynamic calculation 
of torque. These methods would work fine for 

providing simple pump control function, such as 

15 shutting down the pump when it is pumped off. 

However, these methods would not be suitable for real 
time closed- loop pump control, such as rod load 
limiting, that requires a high bandwidth feedback 
signal . 

20 [0009] Past work involving the analysis of rod pump 

systems can be divided into two categories. One such 
category involves predicting the performance of a rod 
pump unit by calculating surface load from known 
surface position and assumed pump load. An example of 

25 this method for deriving the surface dynamometer card 
from the downhole dynamometer card is disclosed in an 
article entitled "Predicting the Behavior of Sucker- 
Rod Pumping Systems", by S. G. Gibbs, in JPT, July 
1963, pages 769-78, Trans, AIME 228. This uses a 

30 multisection model of the rod string to simulate the 
pump operation. 

[0010] The other category deals with the diagnosis 

of existing pumping installations by determining 
actual pump conditions from measured surface 
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conditions. United States Patent No. 3,343,409 

discloses a method for estimating the downhole 
dynamometer card from the surface dynamometer card 
using frequency based Fourier analysis. However, this 
method requires a large number of coefficients to 
accurately model the high frequency components that 
produce the corners of the dynamometer card. In 
addition, the method relies on external sensors for 
polished rod load and position. 

[0011] The average output flow rate of a sucker rod 
pump is a function of the downhole pump stroke and the 
average speed of the pump. With existing 

technologies, the downhole stroke of the pump is 
dictated by the speed of the pumping unit and the 
given characteristics of the pumping unit geometry and 
the sucker rod stiffness. Significant stretch in the 
sucker rod, particularly for deep wells, reduces the 
amount of surface rod stroke that can be delivered to 
the downhole pump. Additionally, the speed of the 
pumping operation is often limited by the need to 
avoid overstressing the sucker rod and/or the pumping 
unit gearbox. Therefore, output flow rate is 

constrained by the imposed pump stroke and stroking 
rate . 



SUMMARY OF THE INVENTION 
[0012] The disadvantages and limitations of the 

background art discussed above are overcome by the 
present invention. With this invention, there is 
provided a method of continuously determining 
operational parameters of a rod pump used in oil and 
gas production, wherein the rod pump includes a rod 
string carrying a downhole pump, the rod string 
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including a polished rod, and a drive system including 
an AC electrical drive motor having a rotor coupled to 
the rod string through a transmission unit. The 
method comprises the steps of continuously measuring 
the electrical voltages applied to the drive motor to 
produce electrical voltage output signals; 
continuously measuring the electrical currents applied 
to the drive motor to produce electrical current 
output signals; deriving values of instantaneous 
electrical torque from the electrical voltage output 
signals and the electrical current output signals; 
deriving values of instantaneous motor velocity from 
the electrical voltage output signals and the 
electrical current output signals; and using geometry 
of the rod pumping unit and one of the instantaneous 
values to calculate instantaneous values of an 
operating parameter of the rod pump. In one 

embodiment, the method is used for calculating rod 
load and/or rod position of a rod pump. The method 
also provides calculations of other pump parameters 
such as gearbox torque and pump stroke that are useful 
in protecting the pumping mechanism and diagnosing 
pump problems. 

[0013] The invention provides a method of deriving 
operating parameters, such as rod load and position, 
from the drive motor and pumping unit parameters 
without the need for external instrumentation such as 
down hole sensors, acoustic fluid level sensors, flow 
sensors, etc. The method provides nearly 

instantaneous readings of motor velocity and torque 
which can be used for both monitoring and real-time, 
closed-loop control of the rod pump. In addition, 
American Petroleum Institute specification geometry 
and system identification routines are used to 
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establish parameters used in calculating the 
performance parameters that are used in real time 
closed loop control of the operation of the rod pump, 
obviating the need to create large look-up tables for 
5 parameter values used in calculating performance 
parameters. Simple parameters defining the special 
geometry used in belt driven pumping units are also 
included in the control. 

[0014] In one embodiment, wherein the first and 

10 second operating parameters are instantaneous position 
and load of the polished rod, the method includes the 
steps of using the estimated values of position and 
load for the polished rod to obtain a surface 
dynamometer card for the rod pump, and deriving from 
15 the surface dynamometer card the instantaneous 
position and load of the downhole pump for pump 
control and/or generation of a downhole dynamometer 
card for the pump . 

[0015] The parameter estimator reduces the cost and 
2 0 complexity of rod pumping systems and provides rod 
load measurement accuracy superior to systems using 
sensors such as strain gages and load cells. 
Moreover, this eliminates wires to sensors mounted on 
moving portions of the pump and reliability issues 
25 related to the sensors and their associated wiring. 

[0016] Further in accordance with the invention, 

the parameter estimator produces values of rod pump 
parameters which can be used in optimizing the 
operation of the rod pump. Thus, in accordance with a 
30 further aspect of the invention, there are provided 
several methods of controlling the rod load and/or 
flow rate of a rod pump used in oil and gas production 
and/or preventing damage to the pump assembly, wherein 
the rod pump includes a rod string including a 
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polished rod and a drive system including an AC 
electrical motor having a rotor that is coupled 
through a transmission unit to the rod string for 
reciprocating a downhole pump. 
5 [0017] One method for rod load control uses the 

computed rod load to control the force in the rod and 
thereby prevent damage to the rod string due to 
excessive tension or compression of the sucker rod. 
Increased pump speeds will typically produce large 

10 tensile force excursions on the up stroke and large 
compressive forces on the downstroke. The method 
limits those excursions to preset limits by 
manipulating the pumping speed. A second aspect of 
the method provides for intentionally increasing or 

15 decreasing rod load during certain portions of the 
pump cycle to increase pump stroke and associated 
fluid production. 

[0018] Another method of rod pump control provides 
for the use of a model of the rod string to derive a 

2 0 factor for modulating pump speed that reduces rod peak 
loads, damps rod force excursions, reduces gearbox 
torque loading, increases pump stroke, and improves 
energy efficiency without the need for external rod 
load and position sensors. Several embodiments of 

25 this method use somewhat different models for control 
of the pump. Those models include the use of rod load 
and/or rod position to generate control signals that 
manipulate pump operation. 

[0019] The rod pump control method comprises the 
30 steps of obtaining a measure of the velocity of the 
polished rod in real-time; obtaining a measure of 
polished rod load in real-time; obtaining an estimate 
of the velocity of the pump in real-time; deriving a 
modulating factor from the difference between the 
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velocity of the polished rod and the estimated pump 
velocity; and using the modulating factor to modulate 
motor speed to cause the downhole pump to more closely 
follow the polished rod position without excessive 
5 excursions in rod load. 

[002 0] The invention allows the stroke of the 

downhole pump to be increased without an increase in 
overall average pumping speed. This increases well 
fluid production without increasing overall pumping 

10 speed and enables increased output in wells that are 
running at maximum physical capacity of the pumping 
system. Alternatively, the method can maintain well 
output with decreased overall pumping speed, reduced 
rod stress fluctuation, and improved energy 

15 efficiency. 

[0021] In accordance with a further aspect of the 

invention, there is provided a system for continuously 
determining operating parameters of a rod pump used in 
oil or gas production, the rod pump including a rod 

2 0 string carrying a downhole pump driven by an 
electrical drive motor that is coupled to the rod 
string through a transmission unit. The system 

comprises means for determining the torque and 
velocity inputs to the rod pump, means for using the 

2 5 torque and velocity inputs to calculate one or more 

values representing the performance of the rod pump, 
and means for using parameters related to the geometry 
of the rod pump and at least one of said performance 
values to calculate values of at least one operating 

3 0 parameter of the rod pump. 

[0022] The rod pump control reduces peak rod loads, 

prevents compressive rod forces, and dampens rod load 
oscillations thereby reducing rod fatigue and rod 
failure. In addition, the rod pump control reduces 
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peak pump velocity, resulting in less power lost to 
viscous pump friction, increasing pumping efficiency 
and reducing pump wear. Moreover, internal frictional 
losses in the rod are reduced by damping rod 
oscillations, thereby increasing pumping efficiency. 

DESCRIPTION OF THE DRAWINGS 
[0023] These and other advantages of the present 
invention are best understood with reference to the 
drawings, in which: 

[0024] FIG. 1 is a simplified representation of a 

rod pump system including a rod pump control system 
that includes a parameter estimator in accordance with 
the present invention; 

[0025] FIG. 2 is a block diagram of the rod pump 
control system of FIG. 1; 

[0026] FIG. 3 is a block diagram of the parameter 

estimator of the rod pump control system for 
calculating values including gearbox torque, polished 
rod load, and rod position using parameters of the 
drive motor and rod pumping unit in accordance with 
the present invention; 

[0027] FIG. 4 is a block diagram of a process for 

obtaining an estimate of rotary weight torque for the 
process of FIG. 3 ; 

[0028] FIG. 5 is a block diagram of a process for 
obtaining an estimate of total reflected inertia for 
the process of FIG. 3; 

[0029] FIG. 6 is a block diagram of a process for 

obtaining an estimate of rod load for the process of 
FIG. 3; 

[0030] FIG. 7 is a block diagram of a process for 

selecting rod stroke regions that have torque factors 
of sufficient magnitude to produce accurate 
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measurement of rod load as the rod load windowing of 
FIG. 3; 

[0031] FIG. 8 is a process flow chart for 

calculating polished rod load and rod position for the 
pump system of FIG. 1, in accordance with the present 
invention; 

[0032] FIG. 9 illustrates a simulated surface and 

downhole dynamometer card for a conventional beam pump 
as well as a downhole dynamometer card generated by 
the wave method of computation; 

[0033] FIG. 10 illustrates simulated surface and 

downhole dynamometer cards for a conventional beam 
pump from a commercially available rod pump simulation 
program; 

[0034] FIG. 11 illustrates a measured surface 

dynamometer card for a belt driven pump as well as a 
downhole dynamometer card generated by the wave method 
of computation; 

[0035] FIG. 12 illustrates simulated surface and 

downhole dynamometer cards for a belt driven pump from 
a commercially available rod pump simulation program; 

[0036] FIG. 13 is a block diagram of a system for 

estimation of rod pump surface and downhole 
dynamometer cards for the rod pumping unit in 
accordance with the present invention; 

[0 037] Fig . 14 is a block diagram of a rod load 

control in accordance with the present invention; 

[0038] FIG. 15 is a block diagram of a single 

section simulation model based rod pump control in 
accordance with the present invention; 

[0039] FIG. 16 is a block diagram of a multisection 

simulation model based rod pump control in accordance 
with the present invention; 
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[0040] FIG. 17 is a block diagram of a wave 
equation based rod pump control in accordance with the 
present invention; 

[0041] FIG. 18 is a process flow chart for 

5 producing rod pump control for improved operation of 
the rod pump system of FIG. 1. 

[0042] FIG. 19 is a surface dynamometer card for a 

beam pump running without rod pump control; 
[0043] FIG. 20 is a surface dynamometer card for a 

10 beam pump running with rod pump control; 

[0044] FIG. 21 is a downhole dynamometer card for a 

beam pump running without rod pump control; 
[0045] FIG. 22 is a downhole dynamometer card for a 

beam pump running with rod pump control; 

15 [0046] FIG. 23 is a graph showing pump velocity as 

a function of time for a beam pump running without rod 
pump control ; 

[0047] FIG. 24 is a graph showing pump velocity as 

a function of time for a beam pump running with rod 
2 0 pump control; and 

[0048] FIG. 25 is a block diagram of a processor of 

the rod pump control system of FIG. 2. 

Definitions of technical terms 
25 The following are definitions of some of the 

technical terms used in the detailed description of 
the preferred embodiments. 

Beam Weight (Wb) : The equivalent weight of the beam 
30 that is used to calculate its articulating 

inertia . 

Counterweight Angle (At) : The angle of the crank 
counterweight offset . 
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Counterweight Inertia (Jc) : The effective inertia of 

the counterweight. 
Crank Angle (Ac) The angular position of the 

beam pump crankshaft at the output of the 
5 reduction gearbox with respect to a reference 

point . 

Crank Velocity (Wc) : The change in crank angle as a 

function of time. The time derivative of the 
crank angle. 

10 Downhole Pump Velocity (Vp) : The velocity of the 

downhole pump as determined by the rod 
string/pump simulation algorithm. 
Electrical Torque (Te) : The torque generated at the 

motor shaft as determined from the motor voltages 
15 and currents. 

Excitation Frequency (We) : The fundamental frequency 
of the instantaneous current circulating in the 
drive motor. 

Gearbox Output Torque (Tn) : The torque at the output 
2 0 of the gearbox. 

Motor Inertia (Jm) : The inertia of the motor and 
associated components rotating at the motor 
speed . 

Motor Velocity (Wr) : The feedback velocity of the 
2 5 motor as determined from the motor voltages and 

currents . 

Overall Gear Ratio (Ng) : The gearing reduction 

between the motor output shaft and the crank 
shaft of the pumping unit . The pumping unit gear 
30 ratio. 

Rod Load (Fr) : The load applied to the polished rod 

as determined by the motor torque, pumping unit 
geometry, and pumping system parameters . 

Rod Position (Xr) : The position of the polished rod 
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as determined by the motor position and the 

pumping unit geometry. 
Rod Velocity (Vr) : The velocity of the polished rod 

as determined by the motor velocity and the 
5 pumping unit geometry. 

Rotary Weight Torque (Tr) : The torque component seen 

at the gear box output shaft due to the 

counterweight normal force. 
Torque Command (Tc) : The final torque command to the 
10 drive system controlling the pump motor. 

Torque Factor (Tf) : A factor that, when multiplied by 

the load at the polished rod, gives the torque at 

the crankshaft of the pumping unit reducer. 
Total Reflected Inertia (Jt) : The inertia seen at the 
15 motor shaft consisting of motor inertia and 

associated high speed components and the 

reflected inertias of the counterweight mass and 

beam mass. 

Unbalanced Force (Bu) : The force that would be 

2 0 required to bring the beam of the pumping unit 

to a horizontal position if the unit had no 
counterbalance . 
Velocity Request (Wx) : The pumping unit prime mover 
operator requested run speed. 
25 Velocity Command (Wy) : The pumping unit prime mover 
command velocity. This signal is a conditioned 
version of the operator requested run speed, and 
originates in the drive control software. 

3 0 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0049] Referring to FIG. 1, there is shown a rod 

pump system 20, the operation of which is controlled 
by a rod pump control system and method including a 
parameter estimator in accordance with the present 
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invention. For purposes of illustration, the rod pump 
control system 21 is described with reference to an 
application in a rod pump system 2 0 that includes a 
conventional beam pump. The beam pump has a walking 
5 beam 22 that reciprocates a rod string 24 that 
includes a polished rod portion 25. The rod string 24 
is suspended from the beam for actuating a downhole 
pump 26 that is disposed at the bottom of a well 28. 
However, the rod pump control system and method 

10 provided by the invention are applicable to any system 
that uses an electric motor to reciprocate a rod 
string, including those that drive the rod through 
belt or chain drives. For example, a belt driven 
pumping unit includes a belt that is coupled to a rod 

15 string for reciprocating the rod string vertically 
within a well as the belt is driven by a motor. 
[0050] The walking beam 22, in turn, is actuated by 
the pitman arm 31 which is reciprocated by a crank arm 
3 0 driven by an electric motor 32 that is coupled to 

2 0 the crank arm 3 0 through a gear reduction mechanism, 

such as gearbox 34. The typical motor 32 can be a 
three-phase AC induction motor operable at 460 VAC and 
developing 10-125 horsepower, depending upon the 
capacity and depth of the pump. Other types of motors 
25 such as synchronous motors can be used to drive the 
pumping unit. The gearbox 34 converts motor torque to 
a low speed but high torque output for driving the 
crank arm 30. The crank arm 30 is provided with a 
counterweight 3 6 that serves to balance the rod string 

3 0 24 suspended from the beam 22 in the manner known in 

the art. Counterbalance can also be provided by an 
air cylinder such as those found on air-balanced 
units. Belted pumping units may use a counterweight 
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that run in the opposite direction of the rod stroke 
or an air cylinder for counterbalance. 

[0051] The downhole pump 26 is a reciprocating type 
pump having a plunger 3 8 attached to the end of the 
5 rod string 24 and a pump barrel 40 which is attached 
to the end of tubing in the well 28. The plunger 38 
includes a traveling valve 42 and a standing valve 44 
positioned at the bottom of the barrel 40. On the up 
stroke of the pump, the traveling valve 42 closes and 

10 lifts fluid, such as oil and/or water, above the 
plunger 3 8 to the top of the well and the standing 
valve 44 opens and allows additional fluid from the 
reservoir to flow into the pump barrel 40. On the 
down stroke, the traveling valve 4 2 opens and the 

15 standing valve 44 closes in preparation of the next 
cycle. The operation of the pump 26 is controlled so 
that the fluid level maintained in the pump barrel 4 0 
is sufficient to maintain the lower end of the rod 
string 24 in the fluid over its entire stroke. 

20 [0052] Referring to FIG. 2, which is a simplified 

representation of the rod pump control system 21 
including parameter estimator in accordance with the 
present invention, the parameter estimator determines 
parameters relating to operation of the rod pump from 

25 motor data without the need for external 
instrumentation. In one embodiment, instantaneous 
motor currents and voltages together with pump 
parameters are used in determining rod position and 
load without the need for strain gauges, load cells, 

30 or position sensors as well as determining pump 
pressure and pump flow without the need for additional 
downhole or surface sensors. The rod position and 
load can be used to control the operation of the pump 
26 to optimize the operation of the pump 26. In 
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addition, American Petroleum Institute (API) 
specifications have been used to define the pump 
geometry that allows the use of readily available data 
from pump manufacturers. System identification 
5 routines are used to establish installation dependent 
parameters specific to the particular pump used in 
calculating performance parameters that are used in 
real-time closed loop control of the operation of the 
rod pump, obviating the need to create large look-up 
10 tables for parameter values used in calculating 
performance parameters . 

[0053] The pump control system 21 includes 
transducers, such as current and voltage sensors, to 
sense dynamic variables associated with motor torque 

15 and velocity. As shown in FIG. 2, current sensors 50 
are coupled to a sufficient number of the motor leads 
for the type of motor used. The current sensors 50 
provide voltages proportional to the instantaneous 
stator currents in the motor 32 . Voltage sensors 52 

20 are connected across to a sufficient number of the 
motor windings for the type of motor used and provide 
voltages proportional to the instantaneous voltages 
across the motor windings. The current and voltage 
signals produced by sensors 50 and 52 are supplied to 

2 5 a processor 53 through suitable input /output devices 
54 . The processor 53 further includes a processing 
unit 55 and storage devices 56 which stores programs 
and data files used in calculating operating 
parameters and producing control signals for 

30 controlling the operation of the rod pump system 20. 

This control arrangement provides nearly instantaneous 
readings of motor velocity and torque which can be 
used for both monitoring and real-time, closed-loop 
control of the rod pump. For example, in one 
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embodiment, computations of motor velocity and torque 
used for real-time, closed-loop control are provided 
at the rate of 1000 times per second. 

[0054] Motor currents and voltages are sensed to 
5 determine the instantaneous electric power level drawn 
from the power source by the electric motor operating 
the well pump. As the rod string 24 that drives the 
downhole pump 26 is raised and lowered during each 
cycle, the motor 32 is cyclically loaded. Depending 

10 on the particular pump installation configuration, the 
walking beam 22 is at a known position during maximum 
and minimum motor loads. The timing of these maximums 
and minimums can define the operational pumping 
frequency and, by integration of the motor velocity in 

15 light of the motor to crank gearing, it is possible to 
estimate the phase position of the pump crank at any 
time. By monitoring the variances of the motor 
currents and voltages as a function of pump crank 
angle, the voltage and current variances can be used 

2 0 together with parameters related to pump geometry to 
calculate estimates of rod position Xr and rod load 
Fr . 

[0055] Referring to FIG. 3, there is shown a block 
diagram of a parameter estimator 2 3 of the rod pump 

25 control system 21 for calculating estimates of 
parameters including rod position Xr, rod load Fr, and 
gearbox output torque Tn. In one preferred 

embodiment, the calculation is carried out by the 
processing unit 55 (FIG. 2) under the control of 

30 software routines stored in the storage devices 56. 

Block 62 responds to signals corresponding to 
instantaneous values of motor current and voltage to 
produce a measure of electrical torque Te of the drive 
motor 32. Block 63 responds to the signals 
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corresponding to instantaneous values of motor current 
and voltage to produce an estimate of velocity Wr of 
the drive motor 32. Block 64 calculates rod position 
Xr and torque factor Tf . Block 65 calculates an 
5 estimate of rod load Fr. Block 67 calculates an 
estimate of rotary weight torque Tr. Block 68 

calculates total reflected inertia Jt . Block 69 
produces an output corresponding to acceleration Alpha 
of the drive motor shaft. 

10 [0056] More specifically, blocks 62 and 63 can 

include hardware circuits which convert and calibrate 
the motor current and voltage signals provided by the 
sensors or transducers 50 and 52 (FIG. 2) into current 
and flux signals. The hardware circuits scale and 

15 translate the current and flux signals into an 
internal frame of reference. After scaling and 
translation, the outputs of the voltage and current 
sensors can be digitized by an analog to digital 
converter. Block 62 combines the scaled signals with 

2 0 motor equivalent circuit parameters to produce a 
precise measure of electrical torque Te . Automatic 
identification routings can be used to establish the 
motor equivalent circuit parameters. Block 63 

combines the scaled signals with motor equivalent 

25 circuit parameters to produce a precise measure of 
motor velocity Wr. 

[0057] In one embodiment, the stator flux is 
calculated from motor voltages and currents and the 
electromagnetic torque is directly estimated from the 
30 stator flux and stator current. Three-phase motor 
voltages and currents are converted to dq 

(direct/quadrature) frame signals using three to two 
phase conversion for ease of computation in a manner 
known in the art. Signals in the dq frame can be 
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represented as individual signals or as vectors for 
convenience. Block 62 responds to motor stator 
voltage vector Vs and motor stator current vector Is 
to calculate a measure of electrical torque Te 
produced by the motor. In one embodiment, the 

operations carried out by block 62 for calculating the 
electrical torque estimate are as follows. The stator 
flux vector Fs is obtained from the motor stator 
voltage Vs and motor stator current Is vectors 
according to equation (2) : 

(2) 

Fs= (Vs-Is.Rs) /s 

( 2A) 

Fds= (Vds-Ids.Rs) /s 

x / (2B) 
Fqs= (Vqs-Iqs .Rs) /s 

where Rs is the stator resistance and s (in the 
denominator) is the Laplace operator for 
differentiation. Equation (2A) and <2B) show typical 
examples of the relationship between the vector 
notation for flux Fs, voltage Vs, and current Is and 
actual d axis and q axis signals. 

[0058] In one embodiment, the electrical torque Te 
is estimated directly from the stator flux vector Fs 
obtained from equation (2) and the measured stator 
current vector Is according to equation (3) or its 
equivalent (3A) : 

Te=Ku. (3/2) .P.FsxIs j^) 
Te=Ku . (3/2) .P. (Fds . Iqs-Fqs . Ids) 

where P is the number of motor pole pairs and Ku is a 
unit scale factor to get from MKS units to desired 
units. 

[0059] In one embodiment, rotor velocity Wr is 
obtained from estimates of electrical frequency We and 
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slip frequency Ws . The inputs to block 63 also are 
the stator voltage Vs and stator current Is vectors. 
Block 63 calculates the motor velocity Wr, In one 
embodiment, the operations carried out by block 63 for 
calculating the motor velocity are as follows. A 
rotor flux vector Fr is obtained from the measured 
stator voltage Vs and stator current Is vectors along 
with motor stator resistance Rs, stator inductance Ls, 
magnetizing inductance Lm, leakage inductance SigmaLs, 
and rotor inductance Lr according to equations (4) and 
(5) ; separate d axis and q axis rotor flux 
calculations are shown in equations (5A) and (5B) 
respectively: 



SigmaLs=Ls-Lm^2/Lr (4) 

then, 

Fr= (Lr/Lm) . [Fs-Is . SigmaLs] (5) 
Fdr= (Lr/Lm) . (Fds-SigmaLs . Ids) (5A) 
Fqr= (Lr/Lm) . (Fqs-SigmaLs . Iqs) (5B) 

[0060] The slip frequency Ws can be derived from 
the rotor flux vector Fr, the stator current vector 
Is, magnetizing inductance Lm, rotor inductance Lr, 
and rotor resistance Rr according to equation (6) : 



Ws=Rr. (Lm/Lr) . [Fdr . Iqs-Fqr . Ids] 

Fdr A 2+Fqr A 2 (6) 

[0061] The instantaneous excitation or electrical 
frequency We can be derived from stator flux 
according to equation (7) : 

We = Fds . sFqs-Fqs . sFds 

Fds^2+Fqs"2 (7) 
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[0062] The rotor velocity or motor velocity Wr can 
be derived from the slip frequency Ws and the 
electrical frequency We according to equation (8) : 

Wr= We-Ws (8) 



[0063] The motor velocity Wr is passed through an 
amplifier 70 and divided by the . gain Ng which 

10 represents the overall gear ratio between the motor 
and the pump crankshaft. Consequently, the motor 
velocity Wr that has been obtained from motor voltage 
and current is converted to crank velocity Wc, which 
reflects the overall pumping unit gear ratio, that is 

15 being produced at the output of the gearbox 34. 

[0064] The crank velocity Wc is integrated in block 
71 to obtain a position which, when combined with the 
overall pumping unit gear and a reference position, 
yields the angular position that is the crank angle Ac 

2 0 of the pumping unit gearbox. The reference position 
can be obtained using a magnetic or optical sensing 
device, a cam limit switch, or similar device, to 
define a reference point within the stroke of the 
pumping unit for each cycle of operation. 

25 [0065] Block 64 calculates the rod position Xr and 
the torque factor Tf using the crank angle Ac obtained 
from the crank velocity Wc, and parameters associated 
with beam pump geometry. As is known, the geometry 
of the pumping unit is defined by the American 

30 Petroleum Institute and can be entered directly into 
the control in that format . One source of API 
specifications is API Specification HE, entitled 
"Specification for Pumping Units", seventeenth 
edition, November 1, 1994. Information entered is 
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dependent upon the class of the rod pump and direction 
of rotation. Typical beam pump parameters that are 
used for calculating the rod position Xr include the 
dimensions of the walking beam, crank radius, and 
5 pitman arm as well as the location of the various 
pivot points in the unit. Those pump parameters are 
readily available from pumping unit manufacturers. 
Simple parameters are also included in the control for 
belt type pump mechanisms that are not specified by 
10 the API standard. Automatic identification routings 
are used to establish installation dependent pumping 
unit parameters such as counterbalance inertia and 
frictional terms. 

[0066] Block 67 combines the crank angle Ac with 

15 the counterweight angle At to produce an estimate of 
rotary torque Tr associated with the weight of the 
counterweight. Referring to FIG. 4, in one 

embodiment, the rotary torque Tr is obtained by 
summing the crank angle Ac with the counterweight 

20 angle At using summing block 72. Block 73 obtains the 
sine of the resultant value. The result is passed 
through amplifier 74, the gain -Mu of which is 
selected to correspond to the counterweight moment, 
producing the rotary torque Tr. 

25 [0067] Referring to FIGS. 3 and 5, block 68 

combines torque factor Tf, produced by block 64, beam 
weight Wb, counterweight inertia Jc and motor inertia 
Jm to produce total reflected inertia Jt . With 
reference to FIG. 5, block 75 obtains the product of 

3 0 torque factor Tf and equivalent beam weight Wb . 

Torque factor Tf is entered twice to square that 
factor. In one embodiment, the gain of amplifier 76 
is 1/G which divides the product by the acceleration 
of gravity such that output of amplifier 76 is the 
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articulating inertia of the pumping unit. The result 
is combined with counterweight inertia Jc at the crank 
in summing block 77 and scaled by amplifier 78, the 
gain l/Ng^2 of which is selected to correspond to the 
5 inverse square of the overall gear ratio. The scaled 
value is combined with motor inertia Jm by summing 
block 79 to obtain the estimate of total reflected 
inertia Jt . 

[0068] Referring to FIGS. 3 and 6, block 65 

10 combines electrical torque Te calculated by block 62, 
rotary weight torque estimate Tr calculated by block 
67, an estimate of total reflected inertia Jt 
calculated by block 68, motor acceleration estimate 
Alpha, produced by block 69, torque factor Tf from 

15 block 67, static friction Sf, crank velocity Wc from 
amplifier 70, viscous friction factor Bf, and 
unbalanced force Bu to produce the rod load estimate 
Fr. With reference to FIG. 6, block 80 obtains the 
product of motor acceleration estimate Alpha and an 

20 estimate of total reflected inertia Jt at the motor, 
the result of which is subtracted from the electrical 
torque Te by summing block 81. This difference is 
scaled by a factor corresponding to the gearbox ratio 
Ng, using amplifier 82. The difference between 

25 electrical torque as modified by motor acceleration 
and load inertia factor, provided by amplifier 82 
minus the static torque Sf, provided by setup testing, 
and viscous torque, provided by multiplying in 
amplifier 84 the crank speed Wc by a viscous friction 

3 0 factor Bf determined during setup, and rotary weight 
torque Tr is divided by the torque factor Tf in block 
85, and the result is summed with unbalanced force Bu 
in summing block 86 to produce rod load estimate Fr. 
A rod load update enable switch 87 and memory element 



1010962 



23- 



88 are used to hold the prior value of rod load at 
around the points where the torque factor Tf goes to 
zero as determined by the rod load update enable 
output of block 60 detailed in FIG 7. 
5 [0069] Referring to FIGS. 3 and 7, block 60 

compares the rod position to the positions defined by 
the ends of the rod stroke Sr determined during setup. 
If the rod position is within the deadzone Dz of 
either end of the stroke, the rod load update enable 

10 output is off and block 65 (FIG. 6) is inhibited from 
updating the rod load Fr. Deadzone Dz is also 
determined during setup. With reference to FIG. 7 , 
relational operator 89 compares the current value of 
rod position Xr to the rod load deadzone value Dz and 

15 outputs a logical true if Xr is greater than Dz. 

Summing block 90 subtracts the rod load deadzone value 
Dz from the rod stroke Sr. Relational operator 91 
compares the current value of rod position Xr to the 
output of summing block 90 and outputs a logical true 

20 if Xr is the lesser value. Logical operator 92 
outputs a logical true only if both relational 
operators 89 and 91 are outputting logical trues. 
[0070] Referring again to FIG. 3, the gearbox 
output torque Tn can be computed from the electrical 

25 torque Te produced by block 62 by using the overall 
gearbox ratio. The value of electrical torque Te 
produced by block 62 is multiplied by a factor related 
to gearbox ratio Ng, using an amplifier 66, to provide 
an estimate of gearbox output torque Tn. 

30 [0071] Block 69 produces an output corresponding to 

acceleration Alpha of the drive motor shaft. One 
method to obtain motor acceleration in Alpha Block 69 
is to differentiate the motor velocity Wr. 
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[0072] Multiplier block 61 produces an output 
corresponding to rod velocity Vr by computing the 
product of torque factor Tf and crank velocity Wc . 
[0073] Referring to FIG. 8, there is shown a 
5 process flow diagram for obtaining estimates of 
polished rod position Xr, polished rod load Fr, and 
gearbox torque Tn derived from the motor current and 
voltage in accordance with the invention. At start- 
up, automatic identification routines are used offline 

10 to estimate various parameters. In one embodiment, 
the automatic identification routines determine 
overall gear ratio Ng and counterweight moment Mu for 
use in further calculations. The overall gear ratio 
is the difference between the motor revolutions and 

15 the crank cycle. The automatic identification 

routines also are used to establish motor equivalent 
circuit parameters as well as installation dependent 
pumping unit parameters, such as static friction 
torque Sf and viscous friction factor Bf . 

20 [0074] Referring also to FIG. 3, after 

initialization, block 93 obtains the instantaneous 
values of motor stator current and motor stator 
voltage from sensors 50 and 52, respectively. As 
described above, blocks 62 and 63 respond to the motor 

25 current and voltage signals from the sensors 50 and 52 
for use in calculating motor torque and velocity. 
Motor stator current and voltage are measured 
continuously allowing the instantaneous values of 
current and voltage to be obtained through the 

3 0 measurement . 

[0075] In block 94, the values of instantaneous 
motor current and motor voltage obtained from the 
measurements are used to derive electrical torque Te . 
In one embodiment, the stator flux is derived from the 
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motor currents and voltages, using equation (2) as 
described above. The electrical torque Te can be 
directly estimated from this stator flux and the motor 
current measured, using equation (3) . 
5 [0076] In block 95, the values of instantaneous 

motor current and motor voltage obtained from the 
measurements are used to derive motor velocity Wr. In 
one embodiment, rotor flux is obtained from the 
measured voltage and current, and stator resistance 

10 and inductance, using equations (4) and (5) as has 
been described. Then, slip frequency is derived from 
the rotor flux, the measured motor current, 
magnetizing inductance, rotor inductance, and rotor 
resistance using equation (6) . An estimate of 

15 electrical frequency is derived from the stator flux 
using equation (7) as described above. Then, motor 
velocity Wr is derived from the slip frequency and the 
electrical frequency using equation (8) as described 
above . 

20 [0077] The motor velocity Wr obtained in block 95 

is used to obtain crank velocity Wc in block 96. In 
one embodiment, the crank velocity is obtained by 
scaling the motor velocity as described above with 
reference to FIG. 3. 

25 [0078] In block 98, the crank angle Ac is obtained 

by integrating the crank velocity Wc obtained in block 
96. A limit switch or similar device may be used to 
determine a reference point within the stroke of the 
pumping unit. The crank velocity Wc is integrated to 

30 get position that combined with the overall pumping 
unit gear ratio and reference position give the crank 
angle Ac. 

[0079] In block 100, rod position Xr is calculated 

using the crank angle Ac together with parameters 
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associated with pumping unit geometry as described 
above with reference to FIG. 3. The use of the system 
parameters with crank angle Ac allows the calculation 
of rod position Xr. 
5 [0080] In block 102, the gearbox torque Tn is 

calculated using the electrical torque Te obtained 
from block 92. The overall gear ratio Ng is also used 
to compute gearbox output torque Tn from motor 
electrical torque Te . 
10 [0081] In block 104, the rotary weight torque Tr, 

is calculated by block 67, the total reflected inertia 
Jt, is calculated by block 68, and motor acceleration 
Alpha is calculated by block 69. 

[0082] In block 106, the combination of the system 

15 parameters and pumping unit geometry with electrical 
torque Te provides the computation of rod load Fr. 
The electrical torque estimate Te is used to obtain 
the rod load estimate Fr. 

[0083] The method of estimating the load and 

20 position of the polished rod at the surface is 
possible without requiring down hole sensors, acoustic 
fluid level sensors, flow sensors, etc. . The values 
of polished rod load and position can be commonly 
plotted in XY format to produce a surface dynamometer 
25 card. The estimation method is a real-time, 

continuously updating method, i.e., it is not 
performed off-line in a batch manner. Moreover, the 
method of estimating a surface dynamometer card for a 
rod pump unit does not employ any load or position 
30 transducers. 

[0084] In accordance with a further aspect of the 

invention, the values of polished rod load and 
position can be used to produce a downhole dynamometer 
card estimate without the need for sensors. Referring 
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to FIG. 13, there is shown a block diagram of a system 
107 for obtaining a downhole dynamometer card without 
requiring down hole sensors, acoustic fluid level 
sensors, flow sensors, etc., using the parameter 
5 estimator 23 described above. The system 107 includes 
a downhole dynamometer card estimator, block 108, that 
uses the polished rod position and the polished rod 
load parameter values obtained by the parameter 
estimator 23 to produce an estimation of the downhole 

10 dynamometer card. Thus, the downhole dynamometer card 
estimation is produced without the need for rod 
position or load transducers. The motor voltage and 
current output signals obtained from measurements of 
the motor voltage and current are used to derive 

15 instantaneous values polished rod load and position 
which are used in producing the estimated downhole 
dynamometer card. 

[0085] The accuracy of the estimation of the 

downhole pump is dependent upon simulating damping 
2 0 forces that are inherent in sucker rod pump systems. 

A viscous damping coefficient is used to model these 
damping forces. 

[0086] More specifically, in one embodiment, an 
estimation of the downhole dynamometer card is 

2 5 obtained using the wave equation to model the force 
trajectory along the rod string in distance and time. 
The wave equation is a linear hyperbolic differential 
equation that describes the longitudinal vibrations of 
a long slender rod. Using the wave equation with 

30 viscous damping, the motion of a sucker rod string can 
be approximated. The wave equation is used only to 
model the rod string and force travelling through it. 
The pump sets the boundary conditions for the wave 
equation at the bottom and the surface prime mover 
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sets the boundary conditions for the wave equation at 
the top. The continuous form of the wave equation 
with constant rod diameter is: 



^ 2 n 2 




v 2 — = — 



+ c — 



(9) 



3x 2 2> t2 



where u is the rod displacement, x is the axial 
distance along the length of the rod, c is the damping 
coefficient, and v is the velocity of force 
propagation in the rods. 



estimating a downhole dynamometer card are disclosed, 
for example, in a paper entitled "An Improved Finite- 
Difference Calculation of Downhole Dynamometer Cards 
for Sucker-Rod Pumps", by T. A. Everitt and J. W 
Jennings, SPE 18189, SPE Production Engineering, Feb. 
1992, pages 121-127. For simplicity, Equation (9) is 
for the case of a constant rod diameter. However, as 
disclosed in the referenced paper of T. A. Everitt and 
J. W Jennings, with modification, this method can also 
account for vairiable rod diameter, including tapered 
rod-strings and rod strings of variable density, e.g., 
steel or fiberglass. Solving the wave equation 

requires only two boundary conditions because only 
steady state solutions are needed. The typical use of 
the wave equation would be to use sampled data of a 
surface dynamometer card from a rod pumping systems to 
do an off-line calculation of the pump downhole 
dynamometer card. In this invention, the wave 

equation is solved on-line for each data point so the 
results can be used in the next sample period for 
control of the pumping system. The two boundary 



[0087] 



Details of the use of the wave equation in 
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conditions are polished rod load Fr and position Xr as 
a function of time. These conditions are produced by 
the parameter estimator 23. 

[0088] The damping coefficient c can be similar to 

5 that presented by T. A. Everitt and J. W Jennings in 
the referenced paper, or that presented in United 
States Patent No. 3,343,409 issued to S . G. Gibbs . 
[0089] The accuracy of the downhole dynamometer 

card estimate can be verified by performing 
10 simulations. One verification procedure that can be 
used is similar to that disclosed in the paper by 
Everitt and Jennings referenced above. 

[0090] Using the multisection simulation disclosed 
in the paper by S. G. Gibbs, referenced above, the 

15 surface dynamometer card load is estimated from a 
given surface position trajectory and pump load and 
position. This method computes new rod position 
estimates in time. Then, using the finite difference 
method disclosed by Everitt and Jennings in the paper 

2 0 referenced above, the downhole dynamometer card is 
estimated from the surface dynamometer card generated 
previously. Then, the estimated downhole dynamometer 
card is plotted with the predicted downhole 
dynamometer card to verify the accuracy of the 

25 estimated downhole dynamometer card. 

[0091] FIG. 9 demonstrates the ability of the wave 

equation method to extract downhole pump operation 
from surface information. An assumed full pump 
condition indicated by reference number 9-1 was used 

30 to simulate surface dynamometer card parameters 9-3. 

Those parameters were used with the wave equation 
method to generate the predicted downhole pump 
dynamometer card 9-2 which closely tracks the 
originally assumed pump dynamometer card. FIG. 10 
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shows the results that were obtained using a 
commercially available simulation program to check the 
results of the multisection simulation and wave 
equation shown in FIG. 9 
5 [0092] The foregoing simulations were conducted for 
a conventional beam type rod pump. However, the 
finite difference method can be used for estimating 
the downhole dynamometer card for other types of rod 
pump units, such as a rod pump unit in which the 

10 driver includes a belt drive. FIG. 11 shows results 
for a pumping unit including a belt that is coupled to 
a rod string for reciprocating the rod string 
vertically within a well as the belt is driven by a 
motor. The graph given by FIG. 11 includes a surface 

15 dynamometer card 11-1 obtained from actual field 
measure data and a predicted downhole dynamometer card 
11-2 obtained from the wave equation method. The 
source data was captured at a constant pump velocity 
of four strokes per minute. 

20 [0093] FIG. 12 illustrates results which are 

similar to those illustrated in FIGS. 11 which were 
obtained using a commercially available simulation 
program. 

[0094] The rod load Fr and/or rod position Xr 
25 parameters obtained using the parameter estimator can 
be used to provide various control functions. By way 
of example, control algorithms can use the rod load, 
rod position, or both to achieve improved pump 
operation. 

30 [0095] Referring to FIG. 14, there is shown a block 

diagram of a system 130 for controlling a pump using 
rod load control. This control algorithm uses the rod 
load Fr, which can be from the estimator in FIG. 3, 
along with maximum upper load and minimum lower load 
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parameters to achieve desirable rod loading. Rod 
loads can be increased in areas of the pump cycle with 
low rod stress to increase pump stroke and associated 
production and/or reduced in areas of the pump cycle 
5 with high rod stress to avoid rod damage. 

[0096] When the torque factor Tf, which can be from 
the estimator in FIG. 3, is positive, the switch 135 
causes the upper portion of the control to be 
selected. Summing block 131 subtracts rod load Fr from 

10 the value Max_Upper__Load, which is determined during 
setup, and outputs the result as Fue . If Fue is 
greater than zero, switch 133 causes it to be 
multiplied by the above upper gain Kau in gain block 
136. If Fue is less than or equal to zero, switch 133 

15 causes it to be multiplied by the below upper gain Kbu 
in gain block 137. 

[0097] Similarly, When the torque factor Tf is zero 

or negative, the switch 135 causes the lower portion 
of the control to be selected. Summing block 132 

2 0 subtracts the value Min_Lower__Load, which is 
determined during setup from rod load Fr, and outputs 
the result as Fie. If Fie is greater than zero, 
switch 134 causes it to be multiplied by the above 
lower gain Kal in gain block 138. If Fie is less than 

25 or equal to zero, switch 134 causes it to be 
multiplied by the below lower gain Kbl in gain block 
138 . 

[0098] Whichever value is calculated is then 
multipled with the absolute value of torque factor Tf 
30 by multiplier block 141. The absolute value of Tf is 
derived by the absolute value block 14 0. The output 
of the multiplier block is added to the velocity 
request Wx by summing block 142 to generate the 
velocity command Wy. 
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[0099] Referring to FIG. 15, there is shown the 

block diagram of a system 110 wherein the rod load Fr 
parameter, which can be from the estimator in FIG. 3, 
can be used along with a one section model of the rod 
5 string based on rod stiffness to provide a control 
function referred to hereinafter as a rod load damping 
control. This control dampens the stress excursions 
in the rod string and causes the downhole pump motion 
to more closely follow the motion of the polished rod 
10 at the surface. Therefore, efficiency and reliability 
of the pump system is increased. 

[0100] Rod load Fr is divided by Rod__St if f ness , 

which is determined during setup, in division block 
111. The result is differentiated by derivative 

15 function block 112 producing a velocity error term. 

If the torque factor Tf, which can be from the 
estimator in FIG. 3, is greater than zero, switch 
block 117 causes the velocity error to be multiplied 
by the gain factor Kup in gain block 113 and then 

20 multiplied by the torque factor Tf in multiplier block 
114. If the torque factor Tf is less than or equal to 
zero, switch block 117 causes the velocity error to be 
multiplied by the gain factor Kdn in gain block 115 
and then multiplied by the torque factor Tf in 

25 multiplier block 116. The result is then added to the 
velocity request Wx by summing block 118 to generate 
the velocity command Wy. 

[0101] Referring to FIG. 16, there is shown the 

block diagram of a system 160, wherein the rod 
30 position Xr and rod load Fr parameters, which can be 
from the estimator in FIG. 3, can be used along with a 
multisection simulation model of the rod string to 
provide a control function referred to hereinafter as 
a simulation model control. 
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[0102] Rod load Fr and rod position Xr are input to 
rod string model block 161. The rod string model 
simulates the rod behavior by dividing the rod string 
into a finite number of elements. Each element has a 
mass and spring constant. The dynamic effects of the 
changing rod load Fr and rod position Xr are 
calculated on each section to determine the velocity 
of the downhole pump. 

[0103] The rod velocity Vr, which can be from the 
estimator in FIG. 3, is subtracted from the pump 
velocity in summing block 162 to determine the 
velocity error term. If the torque factor Tf, which 
can be from the estimator in FIG. 3, is greater than 
zero, switch block 167 causes the velocity error to be 
multiplied by the gain factor Kup in gain block 163 
and then multiplied by the torque factor Tf in 
multiplier block 164. If the torque factor Tf is less 
than or equal to zero, switch block 167 causes the 
velocity error to be multiplied by the gain factor Kdn 
in gain block 165 and then multiplied by the torque 
factor Tf in multiplier block 166. The result is then 
added to the velocity request Wx by summing block 168 
to generate the velocity command Wy. 

[0104] Referring to FIG. 17, there is shown the 
block diagram of a system 170, wherein the rod 
position Xr and rod load Fr parameters, which can be 
from the estimator in FIG. 3, can be used along with a 
wave equation model of the rod string to provide a 
control function referred to hereinafter as a wave 
equation control . 

[0105] The wave equation control is a control 

algorithm capable of damping rod load oscillations, 
reducing rod stress, and increasing pump stroke 
without changing the overall pumping speed, or in the 
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alternative, maintaining the well output with 
decreased overall pumping speed. The wave equation 
control according to the invention increases the pump 
stroke, decreases peaks in rod load and dampens rod 
5 load oscillations. However, average pumping speed is 
not affected. The wave equation control enables 
increased output in wells that are running at maximum 
conventional capability of the pumping system. 
[0106] The wave equation control manipulates motor 

10 velocity to maximize downhole pump stroke. The 
control function provided by the wave equation control 
basically consists of estimating pump velocity state 
by means of a discrete rod string, fluid, and pump 
model. The pump velocity state is then multiplied by 

15 a damping gain and summed with the request velocity. 

This lowers the rod load overshoot through active 
damping while also increasing the downhole pump 
stroke. This results in an increase in output flow 
rate without an increase in overall average pumping 

20 speed which, in turn, increases well output without 
increasing overall pumping speed. This can provide 
increased output in wells that are running at maximum 
capacity. Alternatively, a given well output can be 
maintained with decreased overall pumping speed. 

25 [0107] More specifically, with reference to FIG. 

17, the wave equation control includes a rod string 
model 171 in which the rod load Fr and rod position Xr 
can be from the parameter estimator of FIG. 3. 
[0108] The wave equation control 170 employs a rod 

3 0 string model (i.e., rod string model 171) that 
produces pump velocity Vp and pump position Xp states. 
However in one embodiment, only the pump velocity Vp 
is used in the control function. Although pump 
position Xp is not used for control, pump position can 
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be used to estimate pump stroke Sp . The pump stroke 
information, in turn, can be used to generate flow 
rate information. 

[0109] The rod/pump simulation 171 responds to rod 
5 position Xr and rod load Fr and produces an output 
representative of simulated pump velocity Vp . 
[0110] The rod velocity Vr, which can be from the 
estimator in FIG. 3, is subtracted from the pump 
velocity in summing block 172 to determine the 

10 velocity error term. If the torque factor Tf , which 
can be from the estimator in FIG. 3, is greater than 
zero, switch block 177 causes the velocity error to be 
multiplied by the gain factor Kup in gain block 173 
and then multiplied by the torque factor Tf in 

15 multiplier block 174. If the torque factor Tf is less 
than or equal to zero, switch block 177 causes the 
velocity error to be multiplied by the gain factor Kdn 
in gain block 175 and then multiplied by the torque 
factor Tf in multiplier block 176. The result is then 

20 added to the velocity request Wx by summing block 178 
to generate the velocity command Wy. 

[0111] Referring to FIGS. 18, there is shown a 

process flow diagram for producing simulation model 
control and wave equation control in accordance with 

25 the invention. Block 150 obtains the polished rod 
position Xr. This can be done using the algorithm as 
described above with reference to FIG. 3. 
[0112] Block 152 obtains the polished rod velocity 

Vr. This can be done using the algorithm as described 

30 above with reference to FIG. 3. 

[0113] The downhole pump velocity Vp is obtained in 

block 154. This is obtained using the rod string 
model 161 for the simulation model control or 171 for 
the wave equation control . 
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[0114] Then, the difference of the surface rod 
velocity Vr and the downhole pump velocity Vp is 
obtained in block 156 by subtracting the pump velocity 
from the polished rod velocity, as shown by summing 
5 blocks 162 and 172. 

[0115] The modulating factor is created in block 
158 by applying the damping difference between the 
surface rod velocity and the pump velocity to the 
proportional gain amplifiers selected from 163, 165, 
10 173 and 175 by switch blocks 167 and 177 and then 
multiplying by the torque factor Tf in blocks 164, 
166, 174 and 176. 

[0116] The modulating factor is combined with the 

velocity request Wx by summing blocks 168 and 178 to 
15 produce a command velocity Wy for the drive motor 32. 

The velocity command Wy signal varies as a function of 
the change in rod velocity Vr relative to pump 
velocity Vp . 

[0117] FIGS. 19-24 illustrate results of a beam 

2 0 pump running with and without the simulation model 

control algorithm enabled. Referring initially to 
FIGS. 19 and 20, there is shown a surface dynamometer 
card for a pump running without simulation model 
control and with simulation model control, 
25 respectively. The data is that for a beam pump 
running at seven strokes per minute. 

[0118] As can be seen by comparing the dynamometer 
card in FIG. 19 with the dynamometer card shown in 
FIG. 20, with the simulation model control enabled, 

3 0 the rod stress fluctuation is reduced by lowering the 

peak rod up stroke load while the raising minimum rod 
down stroke load. For example, the dynamometer card 
in FIG 19 shows a peak rod load of about 36,000 pounds 
while the dynamometer card in FIG 2 0 shows a peak rod 
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load of about 33,000 pounds. In addition, the 

dynamometer card in FIG 19 shows a minimum rod load of 
about 13,000 pounds while the dynamometer card in FIG 
20 shows a minimum rod load of about 16,000 pounds. 
5 Rod load oscillation is dampened as can be seen by 
comparing FIG. 19 with FIG. 20. Rod load fluctuation 
of 17,000 (33,000-16,000) pounds with rod pump control 
is 26% less then the 23000 (36,000-13,000) pounds 
without simulation model control. 

10 [0119] FIGS. 21 and 22 show the downhole pump 

dynamometer cards associated with FIG. 19 and FIG. 20 
respectively. As can be seen by comparing the 
dynamometer card in FIG. 21, without simulation model 
control, with the dynamometer card shown in FIG. 22, 

15 with the simulation model control, the pump stroke has 
increased from 255 inches to 282 inches. This 27 inch 
difference translates to an increase in fluid 
production of nearly 11%. 

[0120] Additional advantages of simulation model 

2 0 control can be seen by comparing the graphs in FIGS. 

23 and 24. Graphs in those figures show motor 
velocity Wr, pump velocity Vp, and rod velocity Vr. 
Without simulation model control, FIG. 23, the pump 
velocity reaches peak values that are nearly twice 
25 that of the polished rod and considerable time is 
spent dwelling at zero velocity. When simulation 
model control is enabled, FIG. 24, the peak pump 
velocity more nearly tracks polished rod velocity and 
no time is wasted dwelling at zero velocity. This 

3 0 provides for increased pump stroke without the need 

for high pump peak speeds. 

[0121] In this example, pump stroke is increased 

approximately 11% with no overall change in average 
pumping unit speed. In addition, peak rod load is 
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reduced, minimum rod load is increased, rod load 
oscillation is dampened, and peak pump velocity is 
reduced. 

[0122] Referring to FIG. 25, in one preferred 

5 embodiment, the system provided by the present 
invention, is software based and is capable of being 
executed in a processor 53 shown in block diagram form 
in FIG. 25. In one embodiment, the computer system 
includes input devices 181, such as current and 

10 voltage sensors connected to analog to digital 
converters, output devices 182, such as, a variable 
frequency drive, and a processing unit 55 having 
associated random access memory (RAM) and read-only 
memory (ROM) . In one embodiment, the storage devices 

15 include a database 185 and software programs and files 
which are used in carrying out simulations of circuits 
and/or systems in accordance with the invention. The 
programs and files of the computer system include an 
operating system 186, the parameter estimation engine 

2 0 187, and a control method 188 such as the simulation 
model control engine, rod load control engine, rod 
load damping control engine or wave equation control 
engine, for example. The programs and files of the 
computer system can also include or provide storage 

2 5 for data. The processor is connected through suitable 

input/output interfaces and internal peripheral 
interfaces (not shown) to the input devices, the 
output devices, the storage devices, etc., as is 
known . 

3 0 [0123] Although an exemplary embodiment of the 

present invention has been shown and described with 
reference to particular embodiments and applications 
thereof, it will be apparent to those having ordinary 
skill in the art that a number of changes, 
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modifications, or alterations to the invention as 
described herein may be made, none of which depart 
from the spirit or scope of the present invention. 
All such changes, modifications, and alterations 
5 should therefore be seen as being within the scope of 
the present invention. 



1010962 



40- 



